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Abstract 

We study the interaction of an atom with a quantum guided field in a weakly driven fiber-Bragg- 
grating (FBG) cavity. We present an effective Hamiltonian and derive the density-matrix equations 
for the combined atom-cavity system. We calculate the mean photon number, the second-order 
photon correlation function, and the atomic excited-state population. We show that, due to the 
confinement of the guided cavity field in the fiber cross-section plane and in the space between the 
FBG mirrors, the presence of the atom in the FBG cavity can significantly affect the mean photon 
number and the photon statistics even though the cavity finesse is moderate, the cavity is long, 
and the probe field is weak. 
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I. INTRODUCTION 



In cavity quantum electrodynamics, an interesting regime, called the regime of strong 
coupling, occurs when the maximal atom-field dipole coupling strength exceeds the cavity 
field decay rate and the atomic spontaneous emission rate [1]. In this regime, excitations can 
be exchanged coherently between the atom and the field several times before the incoherent 
decay process occurs, the properties of the field can be significantly modified by the presence 
of even a single atom, and the presence of a single photon in the field can saturate the 
response of the atorn. The effects of single atoms on the cavity field in real time {2, 3| 



have been observed 



Gj. It has been reported that the presence of an atom in the cavity. 



which is tuned to the atomic transition and resonantly driven by a laser field, can lead to 
a dramatic drop in the transmitted intensity [4] . It has been demonstrated |7|, l8| that the 
spatial variation of the cavity mode can lead to a confining potential sufficient to trap an 
atom within the cavity mode even for a single quantum of excitation . Cooling [12| of 

single atoms with single photons in a high-Q cavity has also been investigated. 

Recently, it has been proposed to combine the cavity technique with the nanofiber tech- 
nique to obtain a hybrid system, where the interaction is enhanced by the transverse con- 
finement of the field in the fiber cross-section plane as well as the longitudinal confinement of 
the field between the mirrors. It has been shown that the presence of a fiber-Bragg-grating 
(FBG) cavity with a large length (on the order of 10 cm) and a moderate finesse (about 
30) can significantly enhance the group delay of a guided probe field |l3| and substantially 

lere has been 



15 



-|20|. With 



enhance the channeling of emission from an atom into a nanostructure [ij]. T 
a large body of work involving fiber Bragg gratings over the past two decades 
careful control of the grating writing process and appropriate choice of glass material, a 
FBG resonator can have a finesse of well over 1000 and a linewidth of a few MHz [20]. It 
is worth mentioning that several methods for trapping and guiding neutral atoms outside 
a fiber have been proposed and studied j2ll-l27]. A trapping method based on the use of 
two (red- and blue-detuned) light beams has been studied for large-radius fibers [21| and 



nanofibers 



23| and has recently been experimentally realized |27l |. 



In this paper, we study the interaction of an atom with a quantum guided field in a 
weakly driven FBG cavity. We show that, due to the confinement of the guided cavity field 
in the fiber cross-section plane and in the space between the FBG mirrors, the presence 
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of the atom in the FBG cavity can significantly affect the mean photon number as well as 
the photon statistics even though the cavity finesse is moderate, the cavity is long, and the 
probe field is weak. 

The paper is organized as follows. In Sec. [TTlwe describe the model. In Sec. Illll we derive 
the density- matrix equations for the combined atom-field system. In Sec. [IV]we present the 
results of numerical calculations for the mean photon number, the second-order correlation 
function, and the atomic excited-state population. Our conclusions are given in Sec. |Vl 



II. MODEL 



We consider a two-level atom in the vicinity of a nanofiber with two FBG mirrors (see 
Fig. [T]). The field in the guided modes of the nanofiber is refiected back and forth between 
the FBG mirrors. The nanofiber has a cylindrical silica core of radius a and of refractive 
index ni = 1.45 and an infinite vacuum clad of refractive index n2 = 1. In view of the very 
low losses of silica in the wavelength range of interest, we neglect material absorption. We 
use the cylindrical coordinates (r, if, z), with z being the axis of the fiber. 

In the presence of the fiber, the electromagnetic field can be decomposed into guided and 



radiation modes [28|. In order to describe the field in a quantum mechanical formalism, we 
follow the continuum field quantization procedures presented in {2^. First, we temporally 
neglect the presence of the FBG mirrors. Regarding the guided modes, we assume that the 
single-mode condition is satisfied for a finite bandwidth around the atomic transition 
frequency ojq. We label each fundamental guided mode HEn with a frequency u in this 
bandwidth by an index /i = {u,f,l), where / = +, — denotes the forward or backward 
propagation direction and / = +, — denotes the counterclockwise or clockwise rotation of 



polarization. The quantum expression for the electric positive-frequency component Eg^^ of 
the field in the fiber guided modes is [soj 



Here e*^^-* = e*^^^ (r, (p) is the profile function of the guided mode /i in the classical problem, 
is the corresponding photon annihilation operator, = du is the summation 

over the guided modes, (3 is the longitudinal propagation constant, and /3' is the derivative 



of (3 with respect to u. The constant (3 is determined by the fiber eigenvalue equation 



28|. 



The operators and aj^ satisfy the continuous-mode bosonic commutation rules [a ,,, a | ,] = 



6{uj — uj')5ffiSii'. The explicit expression for the mode function e'^'^^ is given in Refs. 
According to Ref. 



28 



3Q|. 



30[, the rate of spontaneous emission into guided modes is given by 
^^0 V- i_, . r^o,/,0|2 



7gyd 



2eohvg 



(2) 



Here d is the matrix element of the electric dipole moment of the atom and Vg = 1/ ^'{ujq) 
is the group velocity of the guided field. 
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FIG. 1: (Color online) An atom in the vicinity of a nanofiber with two fiber-Bragg-grating mirrors 
driven by a weak guided probe light field. 

Regarding the radiation modes, the longitudinal propagation constant /3 for each fre- 
quency uj can vary continuously, from —k to A;, with k = uo/c being the wave number. We 
label each radiation mode by an index v = {u, f3,m,l), where m = 0, ±1, ±2, . . . is the 
mode order and / = +, — is the mode polarization. The quantum expression for the electric 
positive-frequency component E^j^J of the field in the radiation modes is 



Here e'^'^^ = e*^'^^ (r, Lp) is the profile function of the radiation mode v in the classical problem, 
ay is the corresponding photon annihilation operator, and = duj /^^ dj3 is the 

summation over the radiation modes. The operators ay and aj, satisfy the continuous-mode 



bosonic commutation rules [a,^,a|^,] = 5{uj — uj') 5 
for the mode function e*^'^^ is given in Refs. 
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d — P')5mm'^u'- The exphcit expression 



30|. According to Ref. 30|, the rate of 



spontaneous emission into radiation modes is given by 



7-d = T^E rrf/3|d-e(-'^'-'')(,,^)|2. (4) 

^^0^ ml J-^^ 



Next, we take into account the effect of the FBG mirrors on the mode functions. We 
assume that the two FBG mirrors are identical, having the same complex reflection and 
transmission coefficients R and T, respectively, for the guided modes in a broad bandwidth 
around the atomic transition frequency Uq. In general, we have |-Rp + |Tp < 1, where the 
equality (inequality) occurs for lossless (lossy) gratings. Without loss of essential physics, 
we assume that the gratings are lossless, that is, |-Rp + |Tp = 1. Let the mirrors be located 
at the positions z = ±L/2 along the fiber, where L is the distance between the mirrors. 
The guided modes are modified by the presence of the mirrors. The mode functions of the 
cavity-modified guided modes are obtained, as usual in the Fabry-Perot theory, by summing 



the geometric series resulting from the multiple reflections by the mirrors 3ll-l33|. Inside 



the cavity, the mode functions of the cavity-modified guided modes are given by 

and, hence, the electric positive-frequency component of the field in the cavity- modified 
guided modes is 

= ' E a„e(")e'«/*«-'-"'. (6) 

We assume that the FBG mirrors do not reflect the radiation modes. This approximation 
is reasonable when the distance L between the FBG mirrors is large as compared to the 
fiber radius a and to the wavelength Aq = 27r//co, with ko = ujq/c being the wave number of 
the atomic transition. In the framework of this approximation, the mode functions of the 
radiation modes are not affected by the presence of the FBG mirrors. 

We drive the FBG cavity by a classical probe light field propagating along the fiber in 
a guided mode fip = {up,fp,lp). Let Pin be the incident power. The transmitted power is 
given by 

nut ^.n _ + 4|^|2 sin2 0(^^) " I ') 

Here Q{uj) = (3{u)L + (pR is the phase shift caused by a single cavity crossing and a single 
reflection, with (pn being the phase of the reflection coefficient P, that is, P = \R\e'^'^^. We 
assume that the probe field frequency Up is close to a resonant cavity frequency which 
is determined by the equation Q{(jJ^ = mvr, with m being an integer number. To the first 
order in Up — Uc, we have the expansion Q{ujp) = mn + {L/vg){up — Uc). Hence, we find 

P -P ^^-'^^l'^' (8) 

-Tout -^m/^ I DI2N2 , MlDli/ T /„. \'2/. . , . \2 ' \°/ 



In the framework of the input-output formulation for optical cavities, the evolution of 
;he photon operator a for the field in an empty two-sided cavity is governed by the equation 



34| 



d = -iooctt - y^K/2ai^ + Kj^h^^. (9) 

Here and 6in are the input photon operators for the left and right sides of the cavity, 
respectively, and k is the cavity damping coefficient, which is assumed to be the same for the 
two sides. When the input field a-^^ is an external coherent monochromatic probe field, with 
the frequency cUp and the power Pin = 'hbj^lfl'-^^ai^ , and the input field is in the vacuum 



state, the mean number n = (a'^a) of photons in the cavity is given by 34 1 

Here = Up — ujc is the detuning of the probe field from the cavity resonance and 



is the cavity pumping rate. 

It is clear that the interaction of the quantum cavity field with the external classical 
probe field can be described by the Hamiltonian 

Hp = -ih{r]e''^^*a - H.c). (12) 

Note that the power of the transmitted field is related to the mean intracavity photon 
number n via the formula 34 1 



-Pout = ^hupKn. (13) 
When we insert Eq. ( 1T0|) into Eq. ([T3|) and compare the result with Eq. ([8]), we find 

and 



(1- 






2\R\Lhup 



The external probe field excites the guided cavity mode /ic = (ucfcjc), where fc = fp 
and Ic = Ip- In the single-mode regime, the electric positive-frequency component of the 
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field in the excited guided cavity mode is given by the expression 



E(+) = l./!^a[e(--^-^^)e*^^'^^^ 
2 V gqL 

Hence we find that the interaction between the atom and the quantum guided cavity field 
in the dipole and rotating-wave approximations is described by the Hamiltonian 

Haf = -in{Gaa^ - H.c), (17) 

where a and a"^ are the downward and upward operators for the atomic transitions, respec- 
tively, and 



G = ;j^d 



2 V eohL 



g(<^c,/c,ic)gj/c/3c2 



^^lcV (18) 



is the coupling coefficient. Here r, ip, and z are the cylindrical coordinates of the position 
of the atom. 

In general, the atomic dipole vector d of a realistic atom is a complex vector. We use the 
notation Vq = Vz and V±i = t{Vx ± ^K/)/a/2 for the spherical components of an arbitrary 
vector V. We assume that only one spherical component dq = doi the dipole vector d, where 
q = —1, 0, or 1, is nonzero. Then, we have d • el^'^'f'^^ = [—lYde"^q'''\ Due to the properties 
of the mode profile functions [28], we have e'^^^^'^\r,Lp) = f^~^'^e^'^'''^^'^~'^^e^^^i{r,(f)\. Here 
we have introduced the notation leg"^^! = |ei'^'^'^''| and \e^±}\ = {\ei'^'^'^^\ =F |e^'^'^''|)/v^- 
Furthermore, we assume that the probe field and hence the cavity field are counterclockwise 
circularly polarized, that is, Ip = Ic = +, and that the atomic transition is o"+ polarized, 
that is, q = I. Then, we have G = e-^("^+^)^/2(^M2/g^;^^)i/2 cos {(3^z + rmi/ 2). We 

can remove the phase factor e~*('"+^)'^/^ by performing the transformation ae"*^'"^^-''^/^ = d 
for the photon operators or the transformation cre**^™^^-''^/^ = a for the atomic operators. 
Therefore, we can take the following expression for the atom- field coupling coefficient: 



G'= J^leL'i^l cos(/3,z + m7r/2). (19) 



Note that Q = 2G is sometimes called the vacuum Rabi frequency. 
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In the vicinity of the fiber surface, the atom experiences the effect of the van der Waals 
potential on the internal state energy and on the center-of-mass motion. Let Vg and Ve be the 
van der Waals potentials for the ground state \g) and the excited state |e), respectively. In 
the presence of the fiber, the atomic transition frequency is shifted from the bare frequency 

and is given by 

Uair) = Uo + Veg{r)/n, (20) 

where V^g = Ve — Vg. The Hamiltonians of the atom and the guided cavity field in the case 
of no coupling between them are given by 

Ha = + -ncOaO-z H (21) 

2M 2 2 ^ ^ 

and 

Hp = huca^a, (22) 

respectively. Here p and M are the momentum and the mass of the atom, respectively. When 
we sum up the Hamiltonians f|T2|) . f|T7|) . fl2Tl) . and fl22l) . we obtain the total Hamiltonian 
H = Ha + Hp + Haf + Hp for the combined atom-field system. We decompose H into two 
parts. Ho = hujp{a^a + <Jz/2) and Hi = H — Hq. In the interaction picture, the combined 
atom- field system is described by the Hamiltonian ifjnt = e^^^^/^if/e"*^''*/^, whose explicit 
expression is 

1 

-f^int = ^ - -^l^aOz - hAca^a - ihG{aa^ - aV) 

- ih7]{a - at) + (23) 

Here Aa = Up — Ua is the detuning of the probe field frequency Up from the surface-shifted 
atomic transition frequency Ua- The Hamiltonian ( 123|) is almost same as that for a weakly 
driven microcavity [12]. The difference between the two models is that the presence of 
the van der Waals potential in the case of the nanofiber-based cavity affects the transition 
frequency and the center-of-mass motion of the atom. 

For the treatments in the remaining part of this paper, we assume that the translational 
motion of the atom can be neglected. This situation can be realized in the cases where 
the atom is trapped, the atom is slow (the atomic ensemble is kept at a low temperature), 
or the two- level atom is replaced by a heavy particle like a quantum dot. We note that 
trapping of neutral cesium atoms in a one-dimensional optical lattice above the surface of a 
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nanofiber has been realized experimentally [27| . With the above assumption, we can neglect 
the kinetic energy and the potential in the Hamiltonian fl23l) . Then, we obtain the following 
effective Hamiltonian: 

Hes = —-^AaOz — hAcO'a — ihG{aa^ — aV) 

-ih7]{a-a^). (24) 

We emphasize that the effect of the van der Waals potential on the atomic transition fre- 
quency Ua and the detuning Ac is kept in the above Hamiltonian. It is clear that the effect 
of the van der Waals potential on the atomic transition frequency is negligible in the region 
of large atom-to-surface distances r — a but significant in the region of small distances r — a. 
In our numerical calculations presented in Sec. IIVI we assume for simplicity that the van 
der Waals potential in the case of the fiber is the same as that in the case of a fiat surface, 
that is, Va = —C^a/{f — clY 1 where a = g,e. 



III. DENSITY-MATRIX EQUATIONS 

Let p be the density operator for the combined atom-field system. In the presence of the 
atomic decay and the cavity damping, the time evolution of p is governed by the master 
equation 

p = -[p, Hcs] - -{cr'^crp - 2a po^ + pcrV) 

{a) ap — 2apa) + pa) a). (25) 

Here 7 = 7gyd + 7rad is the decay rate of the atom in the presence of the nanofiber and in 
the absence of the FBG cavity. 

We use the basis |a, n) formed from the internal states |a) = |e), \g) of the atom and the 
number states \n) of the guided cavity field. For the matrix elements Pa,n;a',n' = («, ^\p\(^'j 
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of the density operator p, we find the equations 



Pe,n\e,n' i^(^{n TT-^ Pe,n;e,n' 

- G{Vn' + lpe,n;g,n'+l + Vn+ lpg^n+l;e,n') 

,n—l;e,n' 

- Vn' + lpe,n;e,n'+l - V U + lpe,n+l;e,n') 

- lPe,n;e,n' " " [{n' + n)pe,n;e,n' 

-2^{n' + l){n + l)p 

e,n+l;e,n'+lj ; 
Pg,n;g,n' i^(^(^n ^) Pg,n;g,n' 

+ G'(v^Pg,„;e,n'-l + VnPe,n-l;g,n') 
+ r]{Vn/pg^n;g,n'~l + Vnpg^n-l;g,n' 

- Vn/~TTpg,n;g,n'+l - VnTTpg,n+l;g,n') 
~l~ 'y Pe,n;e,n' ~^ ^) Pg,n;g,n' 

- 2^{n' + l){n + l)pg 

Pg,n;e,n' ^^aPg,n\e,n' i^ci^Tl fT') pg,n;e,n' 



— G{-\/n' + lpg^n;g,n'+l — V^Pe,n-l;e,n', 
+ T]{Vn/pg^n;e,n'-l + V^Pg ,n—l;e,n' 

- Vn' + lpg,n;e,n'+l " VnTTpg,n+l;e,r. 



7 K 

- 2^{n' + l){n + l)pg 

,n+l;e,?i'+l 

1. (26) 



c^Pg^n-^e^n' \_\^ ~^ ^) Pg,n;e,n' 



The master equation fl2^ is equivalent to the following equation for the mean value 
(O) = Tr {Op) of an arbitrary operator O: 

- I {Oa^a - 2a^Oap + a^aO) . (27) 

In particular, we find the equations 

^(a) = tA,{a) + G{a)-'^{a) + r], 

= ^A^{a) + G{aa^)-^{a), (28) 
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and 



— (a^a) = G{a)a + ao"^) — K{a^a) + ri[{a)) + (a)), 



dt 



dt 



(aV) = -^(aV + acr^) - 7(0-^), 
+ 2G{a^aa^ + aV) ^(aV + aa^) 



(ata-aat)+r/((a)-((Tt)). 



(29) 



Here A = — is the cavity-atom detuning. 

In order to get insight into the effect of the atom on the guided cavity field, we use the 
procedures of Ref. to hnearize Eqs. f l28|) and fl29|) . For this purpose, we assume that 
the probe field is so weak that the excited state is hardly occupied and there is at most one 
photon in the cavity. In this case, we have 12 1 



(acr^) = -(a), 
{a)aaz) = —{a)a). 



(30) 



With the help of the above formulae, we can linearize Eqs. ( 12811 and ( 12911 and solve them in 



the steady-state regime. The results are 



12| 



(a) = --^(zA„-7/2), 



\D\ 



\D 



{a) a + ao'^) 



\D\ 



-2i- 



\D\ 



(31) 



where 



D = G^ + K^/A - AcAa - z(A,7 + A,/€)/2. 



(32) 
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It is interesting to note that {a^a) = {a^){a), (cr'^cr) = {(t'^){(t), and (aV) = {a'^){(T). These 
relations are vahd only in the case of a weakly driven cavity. The explicit expressions for 
the mean photon number Nca_-v = (a^a) and the atomic excited-state population = (cr'^cr) 
can be represented in the forms 

and 



- (6-2 + K^/A - A, + ( A,7 + A,k) 2/4 ' ^^^^ 
respectively. Since the parameters G, 7, and A^ vary in space, the intensity of the 
transmitted light depends on the position of the atom. We note from Eq. f p3|) 
that the atom can significantly affect the guided cavity field if the condition 3> 
max(/t7/4, 1 Ac| I Aa|, I Ac|7/2, 1 Aa|/t/2) is satisfied. It has been shown that, in the case of 
a high-finesse microcavity driven by a weak probe field, under the condition of strong cou- 
pling, the position of the atom can be inferred from the intensity of the transmitted light 
2^6]. 

IV. NUMERICAL RESULTS 

In this section, we present our numerical results. We use a truncated basis to solve 
numerically the density-matrix equations fl26l) in the steady-state regime. From the steady- 
state solution for the density matrix p, we calculate the mean number A^^cav = (o^a) of 
photons in the guided cavity field and the population Pe = {cr^cr) of the atomic excited 
state. We also calculate the second-order photon correlation function gi^v = {a^a^aa) / {a^a)'^, 



which characterizes the photon statistics 35|]. It is known that the photon distribution of a 
coherent state is a Poisson distribution, which gives gcd, = 1. When Qcai, > 1 or gcsi, < 1, 
the photon statistics is said to be super- or sub-Poissonian, respectively. The occurrence of 

(2) 

sub-Poissonian photon statistics, indicated by the inequality gcai < 1, means that the state 



of the field is nonclassical 35 1. 



We plot in Fig. [2] the mean number Ncav of intracavity photons and the population 
Pe of the atomic excited state as functions of the distance r — a from the atom to the 
fiber surface. The FBG mirror reflectivity is |i?|2 = 0.9, which corresponds to the finesse 
F = 7r|i?|/(l — l-Rp) = 30, a moderate value. The FBG cavity length is L = 10 cm. 
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FIG. 2: Mean photon number A'cav and atomic excited-state population Pg as functions of the atom- 
to-surface distance r — a. The sohd and dashed Hnes represent the results of the calculations from 
the exact steady-state solutions to the density-matrix equations (I26p and from the approximate 
solutions ([331) and (pll) . respectively. The fiber radius is a = 200 nm. The FBG cavity length is 
L = 10 cm. The FBG mirror reflectivity is \R\'^ = 0.9. The input probe power is Pin = 1 pW (left 
panel) and 5 pW (right panel). The parameters for the two-level atom correspond to the P'2-liiie 
transition QSij2F = 4M = 4 o 6P3/2P' = 5M' = 5 of atomic cesium, with the wavelength 
Ao = 852 nm and the natural linewidth 70 = 5.25 MHz. The frequencies of the probe field, the 
cavity, and the atom in free space are equal, i.e. Up = ujc = ojq. The order m of the cavity resonance 
mode is an even number. The axial position z of the atom corresponds to an antinode of the cavity 
standing- wave field. The van der Waals coefficients are C-^g = 1.56 kHz ^m^ and Cse = 3.09 kHz 



The frequencies of the probe field, the cavity, and the atom in free space are equal, i.e. 
Up = Uc = (jJq. The input probe power is Pjn = 1 pW (left panel) and 5 (right panel) pW, 
corresponding to the mean number n = 0.04 and 0.2, respectively, of intracavity photons in 
the absence of the atom. The order m of the cavity resonance mode is an even number. For 
the parameters of this figure, we find the atom-cavity coupling coefficient G{r, z)\r=a,z=o = 
5.3370, the cavity damping rate k = 7.O270, and the atomic decay rate 7(r)|r=a = l-737o, 
with 7o = 5.25 MHz being the natural linewidth of the D2-line transition of atomic cesium. 
It is clear that the strong-coupling condition 2\G\ > k, 7 is satisfied when the atom is close 
to the fiber and to an antinode of the cavity standing-wave field. The figure shows that 




13 



Nca.v and Pe substantially depend on the atom-to-surface distance r — a. We observe that 
the radial profile of Ncs,v has a minimum and the radial profile of Pe has a maximum. The 
minimum of A^^cav is closer to the surface than the maximum of Pg is. When r — a is large, 
Nca.v decreases and P^ increases with decreasing r — a. Such behaviors are due to the radial 
dependence of the atom-cavity coupling coefficient G. In contrast, when r — a is small, A'cav 
increases and Pe decreases with decreasing r — a. Such behaviors are due to the effect of the 
van der Waals potential on the radial dependence of the surface-shifted atomic transition 
frequency Ua- We emphasize again that the van der Waals potential is negligible in the 
region of large r — a but significant in the region of small r — a. Figure El^d) shows that, 
although the probe field is very weak (Pin = 5 pW) and the mean photon number in the 
absence of the atom is small {n = 0.2), the atomic excitation Pg can reach quite large values 
(e.g. Pe = 0.2 at r — a = 320 nm). Such substantial excitation is a result of the strong 
coupling between the atom and the guided cavity field. Comparison between the solid and 
dashed lines in Fig. |2] shows that the approximate solutions fl33|) and f l34|) almost coincide 
with the results of the exact calculations from the density-matrix equations ( l26ll when n is 
on the order of or smaller than 0.04, that is, when the power Pin of the probe field is on the 
order of or smaller than 1 pW. When the power Pin of the probe field is on the order of or 
larger than 5 pW, the discrepancy between the exact and approximate solutions becomes 
serious. 

The effect of the atom on the quantum guided cavity field is significant only when the 
mean photon number is small. However, it is not easy to measure the transmitted field when 
the mean intracavity photon number is too small. Therefore, we focus on the case where 
the mean photon number in the absence of the atom is less than one but not too small. 
For this purpose, we choose the input probe power Pin = 10 pW, which corresponds to the 
value n = 0.4 for the mean photon number in the absence of the atom. The FBG cavity 
length is L = 10 cm. We plot in Fig. [3] the mean photon number N^a-v, the second-order 

(2) 

photon correlation function ^fcav, and the atomic excited-state population Pg as functions of 
the distance r — a from the atom to the fiber surface. The frequencies of the probe field, 
the cavity, and the atom in free space are equal, i.e. Up = Uc = ujq. Figures [3]^a) and[3][c) 
show that Ncav and Pg substantially depend on the atom-to-surface distance r — a. Figure 
int^b) shows that gcai > 1, that is, the photon statistics is super- Poissonian. It is clear from 
the figure that, in the region of large r — a, where the van der Waals potential is weak, the 
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FIG. 3: Mean photon number A'^cav; second-order photon correlation function g'cav, and atomic 
excited-state population as functions of the atom-to-surface distance r — a. The input probe 
power is Pjn = 10 pW. The FBG cavity length is L = 10 cm. The frequencies of the probe field, the 
cavity, and the atom in free space are equal, i.e. Wp = Wc = wq- The axial position z of the atom 
corresponds to an antinode of the cavity standing-wave field. Other parameters are as in Fig. [2j 

second-order photon correlation function gca!v increases with decreasing r — a. However, in 
the region of small r — a, where the van der Waals potential is significant, gcai decreases with 

(2) 

decreasing r — a. Thus, the r dependence of gcav is basically opposite to that of Nc^v Note 

(2) 

that in the vicinity around the distance r — a = 70 nm, the correlation function ^fcav can 
become larger than 2, that is, the photon distribution of the guided cavity field can become 
broader than the Boltzmann distribution for thermal states. 

We plot in Fig. H] the mean photon number Now, the second-order photon correlation 

(2) 

function gcnv, and the atomic excited-state population Pf, as functions of the axial position 
z of the atom along the fiber. Figure ID^a) shows that Now has minima and maxima at the 
antinodes and the nodes of the cavity standing- wave field, respectively. It is clear that, in 
the fiber axial direction, N^av follows the spatial oscillations of the coupling coefficient G 
[see Eq. (I33l)]. Figure Ht^b) shows that gcai has maxima at the antinodes of the guided cavity 

(2) 

field. The figure also shows that gca.v > 1, that is, the photon statistics is super-Poissonian, 

(2) 

in a broad region around each antinode, and that gca.v = 1, that is, the photon statistics is 
almost Poissonian, in a smaller region around each node. A careful look at the data reveals 

(2) 

that gcai < 1, that is, the photon statistics is sub- Poissonian, when the position of the atom 
is shghtly deviated from a node. The appearance of such a nonclassical state is related 
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(2) 

FIG. 4: Mean photon number A'^cav; second-order photon correlation function ^cav, and atomic 
excited-state population Pg as functions of the axial position z of the atom along the fiber. The 
distance from the atom to the fiber surface is r — a = 200 nm. Other parameters are as for Fig. [3l 

to the fact that the absorption of a photon of a weak field by a ground-state atom (or a 
weakly excited one) can reduce the photon-number spread. Figure Hl^c) shows that the axial 
position dependence of Pg is complicated [see Eq. flMj) ]. Indeed, Pe has shallow minima at 
the antinodes and deep minima at the nodes of the guided cavity field. 

We plot in Fig. [S] the mean photon number A^cav? the second-order photon correla- 

(2) 

tion function gfcav, and the atomic excited-state population Pg as functions of the detuning 
Ac of the probe field from the cavity resonance in the case where the cavity is at exact 
resonance with the atom in free space, i.e. lOc = ujq- Figure [5](a) shows that the spec- 
trum of Nca-v has a vacuum Rabi splitting. The positions of the two peaks are close to 



Avdw/2± a/G^ + K7/4 + where Avdw = Kg/^ is the frequency shift caused by the 

van der Waals potential. Due to the effect of the van der Waals potential, the two peaks 
are not symmetric in height and position. The peak on the negative side of the detuning 
Ac is lower and farther away from the center than the other peak. Comparison between the 
solid lines (for r — a = 200 nm), the dashed lines (for r — a = 100 nm), and the dotted lines 
(for r — a = 50 nm) shows that the smaller the distance r — a, the stronger the asymmetry 
of the peaks. The reason is that the depth of the van der Waals potential increases with 
decreasing distance r — a. Figure Mjo) shows that gcai > 1, that is, the photon statistics 
is super-Poissonian, in a broad frequency region where the corresponding atomic excitation 
Pe is significant [see Fig. [5](c)]. The dashed and dotted lines in Figs. [5](b) indicate that 
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FIG. 5: Mean photon number A'^cav; second-order photon correlation function ^cav, and atomic 
excited-state population as functions of the detuning Ac of the probe field from the cavity 
resonance. The cavity is at exact resonance with the atom in free space, i.e. Wc = wq- The distance 
from the atom to the fiber surface is r — a = 200 nm (solid lines), 100 nm (dashed lines), and 50 
nm (dotted lines). The axial position z of the atom along the fiber corresponds to an antinode of 
the cavity standing-wave field. Other parameters are as for Fig. [3l 

the correlation function S'cav can become larger than 2, that is, the photon distribution of 
the guided cavity field can become broader than the Boltzmann distribution for thermal 
states. Such a broadening of the photon distribution is a consequence of the emission from 
the excited atom. 



0.20 -, 




r-a [nm] 



FIG. 6: Same as Fig. [3] except for the detuning Ac = 30 MHz. 
According to Fig. |5l the effect of the atom on the guided cavity field depends on the 
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FIG. 7: Same as Fig. |4] except for the detuning Ac = 30 MHz. 

detuning of the probe field. To see the contrast between the case of exact resonance and the 
case of substantial detuning, we plot in Figs. |6]and[7]the results of the numerical calculations 
for the case where Ac = 30 MHz. Figure |6]^a) shows that, when r — a is large, A^cav increases 
with decreasing r — a. In this region, the presence of the atom increases the number of 
intracavity photons. Figure H^a) shows that the number of intracavity photons is largest 
when the atom is positioned at an antinode of the guided cavity field. Figures [6]^b) and[71^b) 
show that, outside the nodes of the cavity standing- wave field, we have QcJv < 1, that is, the 
photon statistics is sub-Poissonian. Comparison between Figs. [3] and [6] and between Figs, 
m and [7] shows that the two cases have very different behaviors. 

0.004 -L 




r-a [nm] 



FIG. 8: Same as Fig. [3] except for the cavity length L = 1 mm. 
According to Eq. ( 133|) , the atom can significantly affect the guided cavity field if the condi- 
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FIG. 9: Same as Fig. |4] except for the cavity length L = 1 mm. 
0.004-1 „ 




Ac [MHz] 

FIG. 10: Same as Fig. [5] except for the cavity length L = 1 mm. 

tion ^ max(K7/4, |Ac||Aa|, \ Ac\'~f/2, |Aa|/t/2) is satisfied. In the case where = 0, 

the above condition reduces to ^ kj/A. The strong-coupling condition 2G ^ k, 7 is 
not required. We plot in Figs. ISHTUl the results of the numerical calculations for the case 
where the FBG cavity length is L = 1 mm. Such a length is two orders smaller than the 
length used for the calculations of the previous figures. For the parameters of Figs. [8]-[T0| 
we find the atom-cavity coupling coefficient G{r, z)\r=a,z=o = 53. 870, the cavity damping 
rate k = 7027o, and the atomic decay rate 7(r)|r=a = l-737o. Since the cavity length L is 
small, the atom-cavity coupling coefficient G and the cavity damping rate k are large. It 
is clear that the strong-coupling condition 2G ^ /t, 7 is not satisfied. The reason is that k 
increases with decreasing L faster than G does. Figures [814T0] show that the atom can still 
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affect significantly tfie guided cavity field even in the overdamped regime. Note that the 
shapes of the curves in Figs. [8] and M are very similar to those in Figs. |3] and HI respectively. 
The substantial difference between the magnitudes of the mean photon number A^^^av in the 
case of Figs. [SHTUl and the case of Figs. [3HS] results from the substantial difference in the 
cavity length L. Meanwhile, the shapes of the curves in Fig. [TUlfa) are very different from 
the shapes of the curves in Fig. [5]^a). Indeed, the curves in Fig. [TOT a) do not show the 
vacuum Rabi splitting. Such a splitting can be observed only under the condition of strong 
coupling. 

V. SUMMARY 

We have studied the effect of an atom on a quantum guided field in a weakly driven FBG 
cavity. We have calculated the mean photon number, the second-order photon correlation 
function, and the atomic excited-state population. We have shown that, due to the confine- 
ment of the guided cavity field in the fiber cross-section plane and in the space between the 
FBG mirrors, the presence of the atom in the FBG cavity can significantly affect the mean 
photon number as well as the photon statistics even though the cavity finesse is moderate, 
the cavity is long, and the probe field is weak. Due to the effect of the van der Waals 
potential, the vacuum Rabi splitting can become asymmetric with respect to the positions 
and heights of the peaks. The photon statistics of the quantum guided cavity field can be 
super- or sub-Poissonian depending on the position of the atom and the detuning of the 
probe field. 
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